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We compare theoretical expectations for the Casimir force with the results of precise measure- 
ments. The force is calculated at finite temperature for multilayered covering of the bodies using the 
Lifshitz theory. We argue that the dielectric function of the metallization has to be directly mea- 
sured to reach the necessary precision in the force calculation. Without knowledge of this function 
one can establish a well defined upper limit on the force using parameters of perfect single-crystal 
materials. The force measured in the torsion pendulum experiment does not contradict to the upper 
limit. Importance of a thin Au/Pd layer in the atomic force microscope experiments is stressed. 
The force measured with the microscope is larger than the upper limit at small separations between 
bodies. The discrepancy is significant and reproduced for both performed measurements. The origin 
of the discrepancy is discussed. The simplest modification of the experiment is proposed allowing 
to make its results more reliable and answer the question if the discrepancy has any relation with 
the existence of a new force. 
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I. INTRODUCTION 

> 

The Casimir force |lj (see g for a review) between closely spaced macroscopic bodies is an effect of quantum 
electrodynamics (QED) and for this reason it could be predicted very accurately. The force acting between nonideal 
bodies can be found using the Lifshitz theory where it depends on optical properties of used materials. Knowledge 
of these properties is the weakest element in the theory restricting the accuracy that can be achieved. Experiments 
measuring the Casimir force are of great importance because they are sensitive to the presence of new fundamental 
forces ^,[1 predicted in many modern theories (see and references therein) . To distinguish a new force from the 
background, we should be able to calculate the Casimir force with a precision better than the experimental one. In the 
series of recent experiments this force has been measured with the torsion pendulum (TP) || in the range of distances 
0.6 — 6 fim and with the atomic force microscope (AFM) [lO|-[T^| in the range 0.1 — 0.9 fim. The corresponding 
precisions were 5% and 1%, respectively. 

For two ideal plates the famous Casimir formula jjj for the force per unit area is 

where a is the distance between plates. In the experiments the force is measured between metallized disc and sphere 



because for two plates it is difficult to keep them parallel. In this case (Q) has to be modified with the proximity force 
theorem (PFT) ]13|| , This theorem allows to evaluate the force by adding the contributions of various distances as if 
they were independent and for plate and sphere it is reduced to 

R+a 

F{a) = 2nR J F pl (x) dx, (2) 

a 

where R is the radius of curvature of the spherical surface. The PFT approximation is good for R 3> a that holds 
true in all the experiments. If we use the Casimir expression (Q) for the force in (||), then the force between plate and 
sphere will be 

n , , it 3 he R 
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Eq. (||) was deduced for ideally conducting bodies at zero temperature and three kinds of corrections have been 
considered to take into account their real properties. The correction due to finite metal conductivity was found jl4 15 
on the base of the free electron model, where the optical properties of a metal were described by the only parameter 
Up which is the plasma frequency. The force including corrections up to the second order [jl6|] is 
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(4) 



For typical plasma frequency u> p ~ 10 16 s 1 and separations a < 1 \im the correction will be more than 10%. 
Correction due to finite temperature has been found |p_ T|| for ideal conductors and the resulting force is given by 



Fj (a) = F° (a) ( 1 + ™f (fl 



(5) 



where £ = kBTa/hc « a(/im)/7.61 for T = 300°_ftT. The function /(£) is expressed via an infinite sum but it can be 
represented approximately as fiO = (£ 3 /27r)C(3) — (£ 4 7r 2 /45) for £ < 1/2. The temperature correction is negligible 
for the AFM experiments Jl^,[l2| since £ is small in the important separation range 0.1 — 0.3 fim and is only a minor 
correction in condition of the TP experiment ^j, where the important separation range was 0.6 — 3 /im. The general 
form of the correction due to surface distortions has been found in |18|. If the bodies are covered by distortions with 
characteristic amplitudes A\ and Ai , then the force up to the second order in the relative amplitudes of the distortions 
has the form 



F c d (a) = F°(a) 
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where the functions /i j (x, y) describe distribution of the distortions on the surfaces and (...) denotes averaging over 
the surface area. Corrections due to surface roughness are very important for the experiment JlO[ . 

At first |9|- |Tl|Jl9| ] the experimental data were treated using these corrections to Eq. (||) , but it was realized soon that 
at least the conductivity correction has to be considered on more reliable basis. In more realistic approach the Lifshitz 
theory Q was used to evaluate the force between bodies 12 2Cj. Similar but technically a little bit different method 
was developed in |2l). In these approaches the force depends on the dielectric function of the bodies at imaginary 
frequencies e {iui). It has to be expressed with the dispersion relation via the imaginary part of the function e (u>) on 
the real axis which can be directly measured. However, in any of the experiments the information on e (lu) was not 
collected and the handbook data were used instead. Such data are good only to make an estimate for the Casimir 
force with the accuracy much worse than the experimental one. The reason is that the dielectric function depends in 
substantial degree on the sample preparation procedure as will be discussed below. Nevertheless, it is possible to find 
p3f a reliable upper limit on the Casimir force using only well defined parameters of perfect crystalline materials. In 
this paper we will discuss in detail this limit and its comparison with the existing experimental data. 

The paper is organized as follows. In Sec. || we give a general expression for the Casimir force between sphere 
and plate made of nonideal materials at nonzero temperature. Then, to treat the experimental data, the expression 
for the force is generalized for the case of layered bodies. The choice of dielectric functions and parameters for the 
used materials is described in Sec. [II. In Sec. |y| we define the boundary values of the optical parameters and find 
the upper limit on the force in conditions of each independent experiment. Possible reasons for discrepancy between 
theory and experiment are discussed in Sec. 0. Our conclusions are given in the last Section. 



II. THEORY 



Let us discuss first a reliable way to evaluate the Casimir force in the experimental configurations. The force per 
unit area between parallel plates arising as a result of electromagnetic fluctuations is generalized by the Lifshitz theory 
Q , where the real material is taken into account by its dielectric function at imaginary frequencies e («£) : 
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Here prime over the sum sign means that n = term is taken with the coefficient 1/2 and 
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(8) 



It is supposed that both bodies were made of identical materials. The function e (iCn) cannot be measured directly 
but can be expressed via imaginary part of the dielectric function e" (cj) on the real axis with the dispersion relation 
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e «) - 1 = ^ J du 
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(9) 



Information on e" (to) can be extracted from the data on reflectivity and absorptivity of electromagnetic waves with 
the frequency ui for a given material. 

Applying PFT to Eq. (@) one can find the force between sphere and plate. The integration in (^) can be done 
analytically and we find 



F(a) = 



kTR, 
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(10) 



Special care needs to treat the first n = term. The formal reason is that becomes zero but the integral over 
p diverges. The physical reason is that this term corresponds to the static limit when for metallic bodies e — > oo. 
This means that any parameter characterizing the dielectric function of a metal cannot appear in the n = term in 
contrast with a dielectric for which it will depend on the static permittivity of the material. In the e — > oo limit the 
functions G\^ become — G% = G2 = 1- The formal problem is overcome by introducing the integration over a new 
variable x — 2p£ n a/c and after that one can take £ n = for the n = term. The resulting contribution of the first 
term in the force corresponds to the classical limit F c i (a) for metals 



F c i (a) 
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where £ (n) is the zeta-function. Note that in this limit the force does not depend on the metal parameters as it 
should be for a static field. 

The bare Casimir force (^|) is reproduced from Eq. ( |l0| ) in the limit e — > 00 and T — > 0. The finite conductivity 
correction also can be derived from (^). To this end one considers the limit of small temperature when the sum 
in (|l0|) can be replaced by the integral and supposes that the dielectric function of the metal covering the bodies is 
described by the free electron plasma model. In this model e (iQ is 
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where uj p is the free electron plasma frequency. Typical value of the frequency ui p ~ 10 s is larger than fluctuation 
frequencies C ~ c/a giving the main contribution in (0). Then one can expand the functions in (|ic| ) in powers 
of the parameter C,/oj p and performing necessary integrations one finds exactly the result (Eh for the conductivity 
corrections 1 . In this way the corrections up to the fourth order were found in recent paper |p4| . The temperature 
correction (I) is also reproduced from ([To]) in the limit of ideal metals e — > 00. In this case the linear in temperature 
correction does not survive since the n = term ([ll]) is exactly canceled by the linear in T contribution from the rest 
terms in the sum. As a result the leading correction behaves only as £ 3 . 

The expression ( |To| ) differs from those used in |lj] and |2(J in two respects. First, in the cited papers the integration 
connected with the PFT was not done analytically that complicated numerical analysis. Second, the zero temperature 



1 The correction is actually connected with finite density of free electrons (finite uj p ) since the metal conductivity is still infinite 
for the plasma model. Nevertheless, we will not change the fixed terminology. 
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limit has been taken. This limit was also considered in pl| , though the PFT integral was evaluated explicitly. It 
seems a reasonable approximation at small separations because the temperature correction in ([5]) is proportional to 
£ 3 and, therefore, is small. However, one should remember that this correction was derived in the limit of ideal 
conductors e — * oo. For a real conductor it will be proportional to £ as expected for difference between sum and 
integral and will be important (for details see p5[ ). We have computed the force according to ([h]) and with the 
integral instead of the sum at the smallest separation a = 100 nm tested in the experiments. For the plasma model 
(H) with uj p = 2 • 10 6 s 1 we have found that the difference between the sum and integral is 2.5 pN for T = 300° K. 



It becomes 4 pN for the Drude dielectric function (see Eq. (18) below) with the damping frequency lu t = 5 • 10 s 



These values exceed the conservative estimate for the experimental errors 2 pN |12| . 

In the AFM experiments an additional Auq^PcIqa layer of 20 nm [[To) or 8 nm 11 lj] thick was on the top of 
Al metallization of the bodies to prevent aluminum oxidation. It has to be included into consideration. This layer 
is transparent for the electromagnetic waves with high frequencies ~ c/a since adsorption, proportional to e" (uj), 
is small. For this reason the layer was ignored in [Jl0|[l2|. However, the force depends on e(iQ for which the low 
frequencies dominate in the dispersion relation (||) because of large e" (uj) and that is why we cannot neglect the 
Au/Pd layer. To take it into account, one has to generalize expression for the force (Q) to the case of layered bodies. 
Suppose that the top layer has the thickness h and its dielectric function is £\. The bottom layer is thick enough to 
be considered as infinite and let its dielectric function be e 2 . The method described in [Q for deriving Eq. (|7|) can 
be easily generalized for layered plates. We have to add only the matching conditions for the Green functions on the 
layers interface. After some algebra the result will look exactly as (0) but with more complex Gi. 2 : 

(«i + s 2 ) (p + si) e^ h ' c + ( Sl - s 2 ) (p - Sl ) e -^ h ' c 
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(si + s 2 ) (p - si) e ^ h / c + (si - s 2 ) (p + si) e -<^ h /° 



G (egfi + gigg) (gig + si) e^ h l c + (ggfi - eis 2 ) fciP - si) e~^ h l c 
(e 2 si +eis 2 ) (eip- s x )e^ h / c + (e 2 si - £is 2 ) (eip + si) e~^ s i' l / c ' 

where si j2 are defined similar to s in (^J). The force between plate and sphere is given by ( JToj ) with the above Gi. 2 . 

To see qualitatively the effect of an additional layer, we found the finite conductivity correction up to the second 
order in this case. Than for the force one has 



FZ(a,h)=F?(a) 



1 — 4K(h) — 1- — I K(h)- 



auip 5 \ auJip 



(14) 



where the function K(h) depends on the plasma frequencies of the layers uji p , w 2p and the thickness of the top layer 
h 

, . LJi p + ^ 2 ptanh (hujip/c) 
u>2p + ujip tanh (nu\ p /c) 

When h = the force will depend only on w 2p and in the case h — > oo on u>i p as it should be. The effect of the 
top layer disappears if the plasma frequencies coincide. The top layer will be negligible if hui\p/c <C 1. For typical 
plasma frequencies ~ 10 16 s _1 it is definitely not the case even for h = 8 nm. The opposite conclusion made in {HJ 
was based on the too small value of oj p for gold as will be discussed below (see also pi]). Eq. (jl4|) is not very good 
approximation and was discussed only for qualitative understanding of the effect. For actual calculations we will use 
the exact equations (|l0|), (|l3|). 

Importance of a thin metallic layer on the body surfaces has been stressed first in p6[ . The general expression 
for the Casimir force between layered bodies has been presented in pl| but was not used their for actual calcula- 
tions. Significant role of the Au/Pd layer in the AFM experiments was indicated in our preprint p3 |, where it was 
demonstrated that the effect far exceeds the experimental errors. This conclusion was supported in ^7|, where the 
expressions ( fl3| ) for Gi. 2 were confirmed using a different method to deduce them. However, the authors were uncer- 
tain on applicability of ([l3]) for thin films with ft, < 25 nm because the spatial dispersion of the dielectric function 
can be important for such films. We discuss this effect in Sec. [y] where we argue that the spatial dispersion can be 
neglected because of very short mean free path for the electrons in thin films. 

III. THE DIELECTRIC FUNCTION 

Now we are able to evaluate the Casimir force in real geometry of the experiments if there is information on the 
dielectric functions of used materials: Au, Al, and Au^^PdaA alloy. Strictly speaking, one has to measure these 
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functions in wide range of wavelengths on the same samples which are used for the force measurement. It was not 
done in all of the experiments and to draw any conclusion from them we have to make some assumptions on the 
dielectric functions. At low frequencies Au and Al are well described by the Drude dielectric function ps| ]: 

s = e' + ie", 

e» = l--=^, e"(oj) = ff r (16) 
y ' oj 2 +oj 2 y ' UJ (uj 2 + uj 2 ) v > 

where uj p is the free electron plasma frequency and uj t is the Drude damping frequency. A simple test for validity of 
the Drude model is behavior of the material resistivity 0] which is defined as 

9 M = Im ( — T — 1 ) = , (17) 

\s (l- s(uj))ujJ £ w; 

where Eq is the free space permittivity. The resistivity is frequency independent within the Drude approximation. 
For high-purity single-crystal samples of Au and Al (entries 2 in Table ||) the frequency behavior of the resistivity 
in the infrared range of wavelengths 3 fim < A < 32 /im is shown in Fig. [j]. The data on the dielectric functions 
were taken from pQ], where the results from many original works are collected. The data for e' (uj) and e" (u>) can be 



fitted with (16) to find the parameters uj p and w r . The points and fitting curves for e" (uj) are shown in the same Fig. 
[j]. Palladium definitely cannot be described by ([l6]) since its resistance significantly changes in the infrared range. 
However, it is known experimentally that amorphous metallic alloys can be described by the Drude approximation 
p9[ . The physical explanation for this is associated with large Drude damping of the compounds like Auo.qPcLqa- 

Of course, at higher frequencies when interband transitions are reached the Drude approximation fails. Nevertheless, 
it is very helpful since low frequencies dominate in the dispersion relation. Extrapolating (^) to all frequencies one 
finds 



w2 



£{K) - l + wfwr (18) 

Let us estimate the relative error inserted in (|l^) due to extrapolation. If ujq is the frequency of the first resonance for 
a given metal, then the contribution in e (iQ of the frequency range loq < uj < oo, where the Drude model does not 
valid, will be (uJ P /C) ■ (w T /wo) for £ > ujq. This contribution one can take as an estimate for the absolute error and, 
therefore, for the relative error one has ~ uj t /ujq. For typical values uj t ~ 10 14 s _1 and ujq > 10 15 s _1 the error can 
be as large as 10% but error in the force is smaller. If we will use (|l^) for the force computation and change uj p by 
5% ( 10% correction to e (i() at all frequencies) then the force is changed less than 2%. Moreover, since the interband 
transitions give a correction to ( |l8| ) which is frequency dependent, it reduces the correction to the force further. Of 
course, we can take the interband transitions into consideration exactly using the handbook data in visible-ultraviolet 
range which are not very sensitive to the purity and defect density as it happens in the infrared range. However, we 



are intended to establish the upper limit on the force using to p in (18) which is definitely larger than any real value 
and for this reason we can neglect the interband transitions. 

Therefore, in all cases of interest we can use Eq. (|l^) to describe the dielectric function of a material on the imaginary 
axis. The question is how we should extract the parameters uj p and uj t from the data. We proceeded as follows. The 



data for the complex refraction index n + in = \fe have been taken from |30|. First, the validity of the Drude 
approximation was checked by calculating the frequency dependence of the material resistivity according to (inn). In 
the investigated cases the resistivity is more or less constant in the wavelength range A > 2 fim (u> < 9.4 • 10 14 s _1 ). 
This range gives the most important contribution to the dispersion relation (||) and, therefore, it is the range where 
we have to extract the Drude parameters. Using lo p found from the high frequency region can happen to be wrong. 
For example, sometime the plasma frequency is estimated using the transition point in the reflectivity dependence 
on frequency. It works, not very good though, for Al but gives considerably smaller value for Au than that found 
from fitting e (uj) in the infrared range. Probably such an estimate was taken in (^[[H)) for Au where very small value 
uj p = 3.6 • 10 15 s _1 was used. 

The optimal fitting procedure is described in p8| . The damping frequency is evaluated first from the ratio 
(1 — s') /s" which depends linearly on frequency in the Drude model 

l-e' uj 

-1^ = 17- (19) 

6 UJ r 

After that uj 2 can be extracted from 1 — e' by linear fit. The results together with the statistical errors are collected in 
Table | for those data in |30| which include the optical behavior of Au and Al in the infrared region. This table clearly 
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demonstrates that the Drude parameters depend significantly on the sample which is used to measure the optical data. 
These samples contained different densities of the defects (such as impurity atoms, vacancies, dislocations, etc.) that 
influence their optical properties. In this sense there are no universal material parameters. Reproducible parameters 
one can get only for high-purity single-crystals. In this connection all the attempts to use the handbook data for the 
Casimir force calculation can be considered only as estimates and cannot claim on high precision. 

Actually in any of the experiment we do not know the Drude parameters even with 10% accuracy. That is because 
the optical properties of evaporated or spattered films which cover the bodies can be quite different from those of 
bulk materials and depend on technological details of film preparation. It is known, for example, that the film density 
is typically 0.7 from that of the bulk material if it was not annealed. For the resistivity of spattered and evaporated 
Au |3^] the value po = 8.2 pfl ■ cm has been reported in contrast with the bulk resistivity 2.25 pfl ■ cm. If a metal is 
evaporated or spattered on a substrate, it has a large number of defects. Relatively thick metallic films (> 100 nm) 
are usually exist in polycrystalline form. Defects will reduce the concentration of free electrons n which defines the 
plasma frequency of the material. They also will increase the damping frequency to T and resistivity since the mean 
free path of electrons will shorten. To minimize these undesirable in practical applications effects the films are usually 
annealed at high temperature. In the experiments it was not reported were the bodies annealed or not but 

one can say definitely that it was not done in the AFM experiments because the polysterene ball cannot exist at the 
annealing temperature. Even more defects present in thin films (< 20 — 30 nm) which are usually amorphous. This 
explains why thin films have very large resistivity in comparison with the bulk material. Entries 1 and 4 for AI in 
Table || correspond to the data for thick film samples. They support our expectations that the plasma frequency for 
films should be smaller and the resistivity larger than those for the bulk material. 

IV. THE UPPER LIMIT 

Though we cannot use the handbook data to evaluate the force, one can constrain it for a given experiment. This 
statement is based on the observation that the force ( |Io|) increases every time when lu p increases or lo t decreases. It 
has simple physical meaning: the force becomes larger when the metal reflectivity increases. For us it is important 
that any technological procedures will reduce u) p and increase uj t for a given material. A perfect single-crystal will 
have the largest plasma frequency and the smallest oj t and these parameters arc well defined. One can use them to 
get the upper limit on the Casimir force. The plasma frequency u> p is defined by the concentration of free electrons 
in the metal n and their effective mass m* 

w p = t/^-> ( 2 °) 

where e is the electron charge. For good metals, which we are concerned, m* is close but larger than the electron 
mass. It will be helpful for what follows to use Eq. (|l7]) and instead of the damping frequency uj t take the static 
resistivity p (0) = po as a parameter. The later can be directly measured for any material. 

A. Torsion pendulum experiment 

In the TP experiment || the quartz lens and plate were covered first with Cu of thickness 0.5 pm and then with 
Au of the same thickness. The Au layer is thick enough to be considered as infinite and Cu will not play any role. 
We will find the upper limit on the electron concentration if suppose that every Au atom produce a free electron with 
the mass m* = m e . Then for the Au plasma frequency one finds uip U = 1.37 • 10 16 s^ 1 . The resistivity for crystalline 
gold is Pq u = 2.25 pfl ■ cm. One can compare these parameters with that given in Table Q to make sure that they 
correspond to the limit values. Substituting these parameters in ( |l8| ) and calculating the force according to ( |l0| ) one 
finds the upper limit on the Casimir force F max (a) in the TP experiment. 

To compare the upper limit on the force with the measured force F exp , it is more convenient to consider the residual 
force defined as[j] 

AF(ai) = F ex P (a % ) - F max ( ai ) , (21) 



1 Note that Lamoreaux || used different definition of the residual force F exp (a*) — F° (oj). 
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where Gtj are the separations for which the force has been measured. Theory and experiment will be in agreement if 
AF will not be positive within the experimental errors. The original experimental data were presented for the lens 
curvature radius R = 11.3 cm and the residual force in this case is shown in Fig. ||a. It clearlyindicates the presence 
of some unexplained force at the smallest separations. However, later the author recognized |2q | that he was working 
with aspheric lens which had the curvature radius R — 12.5 ± 0.3 cm in the place where the force was measured. 
The correction was published in erratum ||. The points for AF(ai) with the corrected R are presented in Fig. ||b. 
This time the prediction obviously does not contradict to the experiment but dealing with the upper limit we cannot 
conclude that there is an agreement, either. 

The question about surface distortions in TP experiment has been raised in |33|]. Surfaces of the bodies have not 
been examined in j^] but roughness of the order of 30 — 40 nm is quite typical for a metallic film on a polished substrate 
and correlates with the substrate roughness. Quartz optics is used for near UV light and its surface has to be polished 
with a precision at least A/10, where A ~ 300 nm is the UV wavelength. It supports the value above which is routinely 
observed with atomic force or tunnel microscope. According to (||) the short-scale stochastic distortions give only 
a few percent correction even for the smallest separation a — 0.6 nm. Large-scale deviations seem potentially more 
dangerous |33|l since the correction can be the first order in A/ a especially if we take into account that the lens was 
not spherical. The radius Tint of the interaction area one can estimate as ri n ± ~ \J Ra ~ 1000 am. Therefore, only 
small area on the lens takes part in the interaction. In this place the lens can be represented as part of the parabolic 
surface 

AR, (22) 



2 (R + AR) 2R 2 R 2 

where AR is the error in the curvature radius. Here the second term describes the error in the plate-lens separation 
and it can be taken as the distortion amplitude. This amplitude is maximal for r — ri nt and for the relative amplitude 
one has an estimate 

d.i^i.a.HT'. (23) 
a 2 R y ' 

This value is rather small and according to (|J) the correction to the force will be less but comparable with the 
experimental errors. Moreover, negligible role of the large-scale distortions is actually an experimental fact. The 
region of the plate and sphere used for the force measurement in Q was varied by tilting the lens with the adjustment 
screws and there was no evidence for any variation of the force depending on the region used for the measurement. 

The first attempt to evaluate the Casimir force was undertaken by the author of TP experiment || who takes 
into account the first finite conductivity correction but used very small uj p for Au. Lamoreaux was the first who 
recognized the necessity of more rigorous approach to the force evaluation [ 20| and importance of thin films on the 
metallic surfaces p6[ . His numerical results were not quite good due to the delicate problem with choice of u> p which 
we discussed above. The matter has been settled in |2l],|22j with the result which coincide with ours. However, our 
statement is that the calculated force represents the upper limit but not the force itself. The reason is that evaporated 
Au film will have smaller plasma frequency than 1.37- 10 16 s^ 1 due to large number of defects in the film. To know the 
exact value of the force, one has to measure the dielectric function of the bodies but not to take it from a handbook. 



B. AFM experiment 

Let us discuss now the upper limit on the Casimir force for the AFM experiment |l0| . The plasma frequency for Al 
can be restricted using ( p0[ ) if one supposes that every atom produces 3 free electrons. It gives uj p l = 2.40- 10 16 s _1 that 
coincide with the largest value in Table |[ The resistivity of perfect crystals is pf = 2.65 an ■ cm. Since we successfully 
predicted the plasma frequencies for the best samples of Au and Al, the same way one can use to estimate uj p for 

Au/ Pd. If each Au atom gives one and Pd atom gives not more than two free electrons, then ujp U ^ Pd = 1.69 • 10 16 s _1 . 

This alloy is used in microelectronics and resistivity of the bulk material is known to be p^ w ^ Pd w 30 nfl ■ cm |54| in 
accordance with the statement that alloys have large Drude damping. These data allow to find the upper limit on 
the force using ( jl0| ) with the functions G1.2 defined in (|l3|). 

Before comparing the upper limit with the measured force we have to discuss a few additional aspects concerning 
the experiment. Real surface of the bodies is always distorted and the distortions are especially important to treat 
the data in p0| . The distortion statistics were analysed with the atomic force microscope |19[ |. The force has to be 
averaged over the distorted surfaces and we use for this the procedure developed in |^9| . The major distortions are 
the large separate crystals situated irregularly on the surfaces with a typical lateral size of 200 nm. The height of the 
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highest and intermediate distortions is about hi = 40 nm and hi = 20 nm, respectively. The homogeneous stochastic 
background of the averaged height ho/2 = 5 nm fills the surface between the major distortions. The character of 
roughness on the plate and on the ball is quite similar. The part of the surface occupied by distortions with the height 
hi, hi, and ho/2 was measured as v\ = 0.11, «2 = 0.25, and v$ = 0.64, respectively. These values are the probabilities 
for the corresponding distortion to appear. The body surface is defined in such a way that averaging over distortions 
gives zero result. Then the averaged force is the sum of local forces for all possible kinds of distortions which face 
each other taken with the corresponding probabilities 

2 

pMst {a) = J- v iVj F{ aij ), (24) 

where ay are the local separations defined in [|Tg|] . For us it will be important that the minimal local separation is 
an — a— 54.8 nm. This procedure seems quite reliable but, of course, large distortions give the feeling of uncertainty. 
The upper limit on the Casimir force has to be averaged with the corresponding roughness parameters according to 

(H). 



2J) - in 

The raw force F m measured in the experiments consists of a few components [n0| 



F m = F c {ai + a ) + F e (ai +a ) + C- (ai + a ) . (25) 

Here a\ is the separation from the voltage applied to the piezo corrected to the cantilever deflection, op is the parameter 
chosen in such a way that a = cii + a is the absolute separation between bodies, the first term in (Eq) is the Casimir 
force, the second term is the electrostatic force corresponding to the measured contact potential 29 mV, the third 
term represents the linearly increasing coupling of the scattered light into the photodiodes (see for details). The 
parameters clq and C were determined at large separations, where F c is represented by the bare force (^). Then 
the Casimir force can be extracted from the raw data with the help of (p5|). Similar way to find the Casimir force 
was used in S. Of course, the separation a has to be defined as the distance between averaged surfaces of Au/Pd 
layers. However, the role of these layers have been underestimated in jl(],[lj|. In [fl0[ the Casimir force was found 
from (|2^) but a was interpreted as the absolute separation between Al surfaces [ ^0|jl9| ]. Effectively the Au/Pd layers 
were changed by Al which has larger oj p and, therefore, the force calculated theoretically was overestimated. 

It was indicated JTo| that the thickness of Au/Pd layer is less than 20 nm, that is why for calculations we use 
the conservative value h — 15 nm . This change makes the force only larger. The experimental points from nT 
(triangles) and theoretical upper limit on the force including the roughness correction (solid line) are shown in Fig 
in the small separations range a < 250 nm. If the top layer is changed by Al, it enlarges the force on 15 pN at the 
smallest separation a = 120 nm. It is clear that the top layer definitely cannot be ignored in the force evaluation. 
Variation of u>p l on 10% gives only 1 pN change in the force because of screening effect of the top layer. The same 

variation in u>p U ' Pd changes the force on 2 pN. The resistivity variation of the Au/Pd layer on 30% gives 1 pN effect. 
At larger separations all the effects become smaller. 

We can see from Fig. || that the upper limit is smaller than the force measured at small separations and the difference 
is significant. This conclusion contradicts to that in |19[| , where good agreement between theory and experiment has 
been reported (dotted line) based on the detailed theoretical analysis. We have already stressed the importance of the 
Au/Pd layer but it is not the only reason of deviations. It comes also from poor behavior of the finite conductivity 
correction used in jl9| at small separations a. The correction was based on a simple interpolating formula |l6| for the 
force between two plates 

- — 

pplates ,j = p plates ( a )(l + H_E_) (26) 

\ 3 aujp J 

which is applicable in wider range of separations than the expansion up to the second order (Q). It was used to 
calculate the conductivity correction between sphere and plate. Applying the proximity force theorem to ((2^) one 
gets 

F(a) = 3F?(a) /§(l + !!-£-) " , (27) 



x 4 V 3x auj z 



1 



where the upper limit was moved to infinity since a«£ This integral can be expressed via the Gauss hypergeometric 
function but in [fL9| it was expanded in the series up to the fourth order 
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F(a)=F c °(a) 



c 72 / c V 152 / c \ 532 / c 

4- 



acjp 5 \ ao7 p / 3 \ aui v J 3 \ aw,, 



(28) 



and the equations (|2J|) and ((28|) were used to get the theoretical prediction for the corrected Casimir force (dotted line 
in Fig. ||). We found that the interpolating curve which we got by numerical intergation of ( p7| ) is very close to the 
exact force evaluated according to (|Io|), (||) with the parameters cj p = 1.88- 10 16 s _1 , ui T = from [[l9| and T — 300°^ 
in all range of distances. The small difference between the curves is the temperature effect which disappears when 
T — ► 0. However, the expansion ( p8| ) works very bad at a < 100 nm. Although the separation in the experiment 
exceeds 100 nm, the local distance in ( p4j ) can be as small as 65 nm, where ( p8| ) is absolutely unacceptable. The same 
is true if one uses the expansion up to the fourth order found in |24|] directly form ([l0]). The dashed line in Fig. || 
shows the force calculated according to (|24|) , ( p7| ) with the parameters above. The divergence of the solid and dashed 
curves is the effect of the top layer and different parameters used for Al. The larger uj^ 1 which we are using for the 
upper limit is partly compensated by the top layer and that is why the dashed curve lies not too far from the solid 



C. Improved AFM experiment 

Very important progress has been achieved in |T^| where controlled metal evaporation and smaller thickness of 
Au/Pd layer (h = 8 nm instead of 20 nm) allowed to reduce the surface roughness to the level when the correction 
to the force becomes practically unimportant. Also the contact potential has been considerably reduced and the 
parameter do defining the absolute separation of the surfaces has been found in independent electrostatic measure- 
ments. Interaction between metallized ball and corrugated plate has been probed in It is not a subject of our 
consideration here. The data for flat plate and sphere in this work are actually in very good agreement with that given 
m H and we will not discuss them specially. The roughness parameters have been reduced to the following pl2| : 
hi = 14 nm, v\ = 0.05; hi — 7 nm, V2 = 0.11; ho/2 = 2 nm, vq — 0.84. Unfortunately, the unjustified assumption 
that the force is insensitive to the presence of Au/Pd layer has been inserted in the procedure of the Casimir force 
extraction from the raw data and the following relation has been used instead of ( ^B|) 



F rn = F c (a + 2h) + F e (a) +Ca. (29) 

For this reason we cannot directly use the data in to compare with the theoretical prediction. It becomes obvious 
if we plot in the same figure the measured force in the experiments [jl^| and |l^| (see Fig. |^). One would expect that 
for thicker Au/Pd layer the force has to be smaller, but the actual relation is opposite and the difference is large. 
Fortunately, it is easy to restore the right data. In |T^ ] ao was found from an independent electrostatic measurement 
and the constant C was determined at large separations when the shift on2h — 16 nm in F c argument was practically 
unimportant. The measured Casimir force F c - m was expressed as 

F c _ m (a + 2h) = F m (a) - F e (a) - Ca, (30) 

but the points in Fig. [| taken from [O] represent the force as a function of true separation. Therefore, the force 
presented in the figure was calculated as 

F c - m (a) = F m {a ~ 2h) - F e (a - 2h) - C{a - 2h). (31) 

The Au/Pd layer certainly cannot be ignored and the right expression for the measured Casimir force must be 

F c _ m (a) = F m (a) - F e (a) - Ca. (32) 

It is obvious that to restore the right data one has to shift the open squares in Fig. | on 2/i = 16 nm to larger 
separations. After this shift a good agreement between two different experiments is reached. To be absolutely sure 
that the right transformation was done we have tried to reproduce the measured force directly from the raw data 
presented in [|l0| and |l^] , where the procedure has been described in details. The data were available only for one scan 
and for this reason our calculations had restricted precision, but it was enough to make a conclusion on reproducibility. 
To check the procedure, we successfully reproduced the force from the raw data in The force found from the 

raw data in [[l2| according to ( |32| ) agrees much better with the shifted points than with the ones presented in Jl2| . 
These detailed explanations are given not only to answer the criticism of our preprint (23| but also because of great 
importance of the conclusion. It is stated in Pa,E3l that the points have to be shifted to smaller separations. Even 
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ignoring the arguments above it is obvious from Fig. ^ that such a shift would give drastic disagreement between two 
experiments made by the same method. 



The upper limit on the Casimir force in conditions of the experiment |12| one can find exactly as was explained 
above. The only difference is the other set of roughness parameters, but in this case the roughness correction is on 
the level of experimental errors. The experimental points from |l2| shifted on 16 nm as was discussed above and the 
corresponding upper limit (solid line) are presented in Fig. ^|. Again we can see that the upper limit is smaller than 
the measured force and deviation increases at smaller separations. Moreover, even if we replace the 8 nm thick top 
layer by Al with the maximal plasma frequency u> p = 2 AO ■ 10 16 s _1 , this disagreement will not disappear as shows the 
dashed line. The residual force defined according ( pl| ) in the experiments jlO| (triangles) and Q| (open squares) is 
shown in Fig. ^. It clearly demonstrates the presence of some unexplained attractive force which is decreasing rapidly 
when the distance between bodies increased. The points from two different experiments are in reasonable agreement 
with each other that means that the residual force is reproducible. The residual force becomes larger if we deviate 
the parameters from their limit values but the agreement between two experiments is not broken. 

V. DISCUSSION 

Let us discuss now possible reasons for disagreement between experiment and theoretical expectations. As was 
mentioned above the main problem is the values of the material parameters which can significantly deviate from their 
handbook values for evaporated or spattered metallic films. The idea of this paper was to find the upper limit on the 
force instead of the force itself. It allowed to use only well defined parameters of perfect single-crystal materials. We 
took the largest values for the plasma frequencies and the smallest ones for the resistivities. Any possible deviation 
from these values will make the force only smaller and disagreement between theory and experiment will be larger. 

Some doubts were raised |27| about the possibility to describe the thin top layer by a dielectric function which 
depends only on frequency. It was stated that the spatial dispersion can be important for thin films because the 
distance traveled by electron during one period of the field can be larger than the film thickness 

— >h, (33) 

where vp is the velocity of the electron on the Fermi surface. This dimensional effect is really exist (see, for example, 
but it is difficult for observation at room temperature. The reason is that for thin metallic films the mean free 
path for electrons is very short (< 100 A) because of large concentration of the defects. Typically the resistivity of 
very thing films is on the level of 100 ^Vt ■ cm. Then for uj p ~ 10 16 from ( |l7| ) one finds ui T ~ 10 15 s -1 . The mean 
free path is estimated as 

I = „ io 1 (34) 

that is smaller than the used Au/Pd film thickness and, therefore, the spatial dispersion can be neglected. That is 
why the standard dependence for the dielectric function e (u>) is widely used in optics of metals up to the film thickness 
in a few nanometers when quantum effects become involved. For the same reason one can neglect the anomalous skin 
effect for evaporated (spattered) films even for thick ones. Extremely high resistivity 2000 ufl ■ cm for 60 nm thick 
Au/Pd film has been reported in p7[ . However, the authors themselves stress that the resistance of the film is, to all 
appearance, dominated by grain boundaries but optical properties of the film are quite usual. This example shows 
once more that the details of the spattering technology cannot be ignored. 

uncertainty was expressed about applicability of the proximity force theorem. At the moment there is no any 
work where the force between sphere and plate was calculated from " the first principles" . There are some heuristic 
approaches Q allowing to calculate nonadditive Casimir force which agree well with the result found by using PFT 
(see additional discussion and references in p7|). The PFT states that the main contribution to the force can be 
found by adding the contributions of various distances as if they were independent and it is applicable to nonadditive 
forces. An example is the electrostatic force which is nonadditive because the surface charge density is nonuniform for 
curved surfaces. One can check that the PFT gives in this case the correct result with the accuracy ~ a/R. For the 
discussed experiments the correction is very small (a/R ~ 0.001). Even if this term appears with a large coefficient, 
say ~ 10, the correction will be only on the level of the experimental errors. 

In the AFM experiments the electrostatic attraction between bodies because of contact potential was carefully taken 
into account. Of course, the aluminum surfaces were partly oxidized and electrons could be trapped in the oxide. 
These charges can be potentially dangerous if the Au/Pd film is not continuous. In this case the trapped charges 
and their images in underlying aluminum will be the source of the dipole field. Then an additional force can arise as 
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a result of dipole-dipole interaction. However, it is difficult to make a reliable estimate for this effect because we do 
not know the concentration of trapped charges and the size of islands in Au/Pd layer or even do discontinuities exist 
at all for the used layer thickness (it depends on details of the covering procedure). In this connection to make the 
experiment absolutely clear, it is preferable to use Au instead of Al metallization because its non-reactive surface has 
strong advantage over Al. It excludes also additional uncertainties connected with Au/Pd layer. One can use as well 
silver or copper but they are not as inert as gold. It is difficult to measure the dielectric function at the wavelengths 
larger than 30 pm but this range gives an important contribution to the dispersion relation. That is why the material 
behavior in this range has to be predictable. One can say definitely that the materials of platinum group cannot be 
used since they are not described by the Drude dielectric function at low frequencies. 

One can speculate that the observed discrepancy is explained by a new Yukawa force mediated by a light scalar 
boson. Then interaction of two atoms is described by the Yukawa potential 

he 

V Y (r) = -aN^— exp(-r/A), (35) 
r 

where a is a dimensionless interaction constant, A is the Compton wavelength of a particle responsible for the 
interaction, and Ni t 2 is the number of nucleons in atoms of the interacting bodies. An additional advantage of Au 
metallization is higher density of the bodies coating. In this case the Yukawa force will be enlarged roughly by the 
factor (pau/ Pai) 2 ~ 50, where pau,ai are the material densities. If the observed discrepancy has relation with the 
Yukawa interaction, the AFM experiment with Au metallization of the bodies will definitely reveal this new force 
even without detailed knowledge of optical properties of the metallization. 



VI. CONCLUSION 



We have analysed the results of recent precise measurements of the Casimir force using the Lifshitz theory to evaluate 
the force. Layered structure of the bodies coating was taken into account in the frame of the Lifshitz approach. It was 
stressed that the force cannot be predicted with necessary accuracy if there is no detailed information on the dielectric 
function of the bodies coating. Fortunately, all the used materials (Au, Al and AuQ^PdoA) are well described in terms 
of Drude parameters uj p and lo t in the infrared range which dominates in the dispersion relation for the dielectric 
function e(iQ). It was noted that one can find the upper limit on the Casimir force that realized for perfect single- 
crystal materials for which electrical and optical properties are well defined. The surface roughness and linear in 
temperature corrections were taken into consideration. It was shown that the upper limit on the Casimir force does 
not contradict to the result of the torsion pendulum experiment 11. The main conclusion of the paper is that the 
upper limit is smaller than the observed force in the AFM experiments and the difference far exceeds experimental 
errors and theoretical uncertainties for small separations between bodies. The simplest modification of the experiment 
is proposed allowing to reveal origin of the discrepancy. 
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FIG. 1. This figure demonstrates validity of the Drude approximation for Al (triangles) and Au (circles) in the infrared 
range. The resistivity does not depend on frequency (left axis) . Solid lines (right axis) show that s" {to) depends on u> according 
to llq) with the parameters given in Table (entries 2). 



FIG. 2. The residual force ( pT| ) as a function of sphere-plate separation a for TP experiment ||. For the original value of 
sphere radius R = 11.3 cm the points of closest approach in (a) demonstrate presence of some unexplained force. With the 
corrected value R = 12.5 cm |2^| the residual force is shown in (b). In this case there is no contradiction between theory and 
experiment. 



FIG. 3. The Casimir force measured in the AFM experiment [jlOj (triangles). The solid line represents the upper limit on the 
force. The dotted line is taken from |lij] where the Au/Pd layer was ignored and expansion (^) used for the finite conductivity 
correction. This expansion fails at small separations and Eq. (^) has to be used instead. The result is represented by the 
dashed line. The difference between the solid and dashed lines is due to Au/Pd layer. 



FIG. 4. Comparison of the results of two AFM experiments. The data from jlOj are marked by solid triangles and the 
data from [l^] presented as the open squares. The obvious contradiction of two experiments made with the same technique is 
connected with unjustified assumption of the transparency of Au/Pd layer in The right data can be restored by simple 

shift all the open squares on 16 nm to larger separations (see explanations in the text) and after that the experiments will 
agree with each other. 



FIG. 5. The data from |l2| shifted on 16 nm (open squares) and the upper limit on the Casimir force (solid line). The 
dashed line represents the force for the case when the 8 nm thick Au/Pd top layer is changed by Al with the maximal plasma 
frequency uj v = 2.4 ■ 10 16 s -1 . 
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FIG. 6. The residual force defined as (£l|) for the two AFM experiments. The data from [jlOj and jl2j are presented as solid 
triangles and open squares, respectively. The figure demonstrates presence of some unexplained force which decreases rapidly 
with separation increase. The residual force is reproduced for both of the experiments. 



Al 


ujp ■ 10 


LU T ■ 10 


po pO. ■ cm 


1* 


1.54 ±0.01 


15.5 ±0.6 


7.39 


2 


2.235 ± 0.001 


12.49 ±0.01 


2.83 


3 


2.43 ± 0.05 


14.4 ±0.7 


2.76 


4* 


1.63 ± 0.03 


18.2 ±0.7 


7.74 


Au 


uj p ■ 10~ io 


LO T ■ 10" 13 


po pQ ■ cm 


1 


1.280 ±0.001 


3.29 ± 0.05 


2.27 


2 


1.3720 ± 0.0006 


4.060 ± 0.002 


2.44 


3 


1.34 ±0.02 


7.08 ±0.18 


4.46 


4 


1.0513 ±0.0007 


6.24 ±0.21 


6.40 



TABLE I. Parameters of the Drude dielectric function ( ]1 6| ) for Al and Au found by fitting the data [gOj of different mea- 
surements in the infrared range (A > 2 pm). The statistical errors of fitting are indicated. The resistivity was calculated via 
ujp and lot according to (|l7|). The stars in the first column mark the data for film samples. 
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